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We study the effects of confinement between two parallel walls on a two dimensional fluid with 
competing interactions which lead to the formation of particle micro-domains at the thermodynamic 
equilibrium (microphases or microseparation) . The possibility to induce structural changes of the 
morphology of the micro-domains is explored, under different confinement conditions and tempera- 
tures. In presence of neutral walls, a switch from stripes of particles to circular clusters (droplets) 
occurs as the temperature decreases, which does not happen in bulk. While the passage from 
droplets to stripes, as the density increases, is a well known phenomenon, the change of the stripes 
into droplets as an effect of temperature is rather unexpected. Depending on the wall separation 
and on the wall-fluid interaction parameters, the stripes can switch from parallel to perpendicular 
to the walls and also a mixed morphology can be stable. 

PACS numbers: Valid PACS appear here 



The formation of spatially modulated phases at the 
thermodynamic equilibrium (microphases or microsepa- 
ration) is a phenomenon occurring in very different sys- 
tems. Examples of 2D microphases can be found in ad- 
sorption films on a solid substrate [l|, Q and in many 
dipolar fluid films 0, 0, S, @, 01 ■ Examples in 3D involve 
colloidal suspensions [p, [9| and amphiphilic systems such 
as block copolymers [101 [T]| . Typical 2D particle do- 
mains include droplets, stripes, strings, rings, bubbles, 
forming either regular or disordered configurations. In 
3D, instead, spherical clusters, array of cylinders and 
lamellae are obtained. Interest in such phenomena is 
very high, also for the number of possible applications in 
nanotechnology and material science, where new devices 
with specific electric, optical and rheological properties 
are designed. For example, as the 2D patterns can be eas- 
ily transferred on different substrates, they can be used as 
masks in nanolitography p^ . Moreover, connections be- 
tween the general mechanisms governing the microphases 
and other branches of science are currently under inves- 
tigation: from vitrification and gelation [l ^ . [l^. [isj . to 
particle aggregation in colloidal solutions to the 

aggregation of proteins attached to biological membranes 
[l7j and to the cluster formation in various disease pro- 
cesses. 

The presence of anisotropic forces often leads to the 
formation of micro-domains. In the last few years it 
has been recognized that also an isotropic interaction can 
lead to microphases. For instance, in colloidal systems, 
which this work is addressed to, the competing interac- 
tions scenario is often invoked to explain the formation of 
micro-domains. This model is based on the use of an ef- 
fective potential, describing the interparticle interactions 
mediated by the environment, which can be a solvent, 
a substrate or an external field. Within this frame, the 
effective potential is made up of a short-range attraction 
(for instance, due to depletion forces) plus a longer-range 
repulsion (for instance, due to partially screened electro- 



static forces): the attraction favouring the condensation 
of particles, while the longer-range repulsion limiting the 
cluster growth, as firstly pointed out in [l^ . The cluster 
morphology depends on the ranges of the competing in- 
teractions and on the density, but also the temperature 
can have a role. It is interesting to notice that other 
forms of the interparticle potential lead to cluster for- 
mation, such as the softened-core models in [l^ and the 
generalized repulsive interactions in [20| . 

While the th eory of simple fluids under confinement is 
well developed [2l|, [24I , the general behaviour of complex 
fluids, such as those undergoing microseparation, has re- 
ceived very limited attention. It is all the same clear that 
confinement is expected to be a powerful tool to modify 
and even induce new pattern morphologies, which might 
not appear at all in the bulk. This has been shown for 
some studies of block copolymers between smooth or pat- 
terned substrates [H, [2J, [2^ [H, |23| . We have no knowl- 
edge, instead, of studies related to the confinement of 
colloidal systems developing microphases, which is the 
situation under investigation in this paper. In particular 
we treat the 2D case, so that this study is relevant, for 
instance, for Langmuir monolayers: an insoluble film of 
particles, entrapped at the interface of two fluids such as 
air and water. Moreover, the use of lateral barriers is a 
typical tool to influence the system, such as modifying 
the surface particle concentration untill the collapse of 
the monolayer. 

In our model the particle-particle potential has an im- 
penetrable core of diameter a, followed by an attractive 
well and a longer-range repulsive tail. The specific form 



we use reads: 

Upp{r) = 

U{r) = 




£acr . r 



if r < cr 

UiRcut) iff! <r<i?e„t (1) 
otherwise 



r being the interparticle distance, Rcut the separation 
at which the potential is cut and shifted to zero. The 
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potential parameters are: Ra — 1 cr, Ca = 1 for the 
short-range attraction and = 2 cr, = 1 for the 
fonger-range repulsion; the cutoff is Rent = 10. Apart 
for a constant term in cclJT]), Upp{r) is the same inter- 
action studied at first in [30j| and then in [l^, H^, [3l[ for 
bulk 2D fiuids. The fluid is confined by smooth parallel 
walls, set at X = ± Lx/2, being the length of the 
simulation box side. Periodic boundary conditions are 
implemented in the y direction. As wall-particle interac- 
tion we assume the following form: 



oo if |x| > 

V{x) — V{x ~ 0) otherwise 



(3) 
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The quantity a is a dimensionless factor, 
Uq = \Upp{r = cr) I, so that the wall-particle 
potential energy for a particle in contact with the wall 
is Vb = V {x = Lx/2 — a/2) ~ a Uq. R^p is the 
range of the wall-particle potential. When a = we 
have two neutral hard walls at distance Lx] a > 
represents repulsive walls. Since we want to study the 
generic behaviour of the system, it is convenient to use 
as control parameter the temperature T. In specific 
cases, the control parameter has a different meaning, like 
the concentration of non adsorbing polymers with the 
depletion interaction. Hereafter, every physical quantity 
is expressed in reduced units: energies and temperatures 
in units of C/o, lengths in units of cr, densities in units 
of CT^, specific heat in units of ks (the Boltzmanu 
constant). Montecarlo simulations are performed in the 
NVT ensemble {N number of particles, V volume, T 
temperature). Specific heat data are obtained through 
the Parallel Tempering (PT) [sl] technique. 

Microphase pattern depends on the particle mean 
density p\ in bulk, the system goes from droplets to 
stripes [29| to bubbles as p increases. In this letter we 
focus our attention mainly on the striped phase. For 
example, at p = 0.4, the bulk fluid with potential ([T][2|) 
develops a grid of parallel stripes, the period of which 
is P = 10.8 ± 1.5. In a previous study [3l| wc have 
shown that the value of the cutoff Rcut lias a negligible 
effect on the period, for instance, an increase of 50% in 
Rcut changes P by less than the uncertainty on the value 
of P. The transition, from the homogeneous phase to the 
striped phase, is signaled by a peak in the specific heat, 
around the temperature T* « 0.58. At lower T, around 
0.4, the particles order on a triangular lattice within the 
stripes. Wc will show that the confined system behaves 
differently for T < 6* or T > 6* with 6 k. 0.4. 

The behaviour of the system at p = 0.4, with neutral 
walls (a = 0), is firstly analyzed for intermediate 
temperatures {6 < T < T*). The general trend of the 
potential energy is to increase with (Fig. [1]) towards 
the bulk value. This is related to the growth of the 
fraction of particles (those in the centre) which do not 
"see" the walls, so that they behave as in the bulk. 




FIG. 1; Potential energy vs wall separation for three densities at 
T = 0.5. The system is confined between two neutral hard walls 
along the x direction. The insets show some snapshots at p = 0.4. 
The vertical line represents the stripe bulk period P. The shaded 
area emphasizes the oscillating profile of U for p = 0.37 and 0.41. 



On the contrary, the particles nearby the walls feel 
a reduced longer-range repulsion, since they interact 
with fewer particles: their energy is therefore more 
negative than in bulk. What is remarkable is that the 
energy profile at fixed T exhibits a number of maxima 
and minima as function of Lx- These extrcma of U 
are connected to different patterns (insets in Fig. [T|) 
stabilized by the walls. For L^. < 12, the stripes 
alignment is perpendicular to the walls. Wc stress that, 
even if Lj. is large enough to fit in one parallel stripe, 
that never happens: the system switches abruptly 
from perpendicular stripes to two parallel stripes. For 
Lx > 15: energy maxima correspond to the most 
frustrated configurations, where two stripes form close 
to the walls and perpendicular ones in the middle; 
energy minima correspond to the unfrustrated parallel 
stripes, occurring for wall separations Lx ~ nP + A {n 
integer, P bulk period, A stripe width). Similar results 
are found also for slightly different densities such as 
p = 0.37 - 0.41 and plotted in Fig. [T] 

Switching on the wall-particle repulsion [a > 0), 
we find a striking dependence of the stripe morphology 
on the parameters of Uwp{x). In this case the particles 
are squeezed in the central region, so that the space 
truly accessible is reduced (right panels in Fig. and 
the system is shifted at an effective density greater 
than the nominal one. We might think to replace the 
repulsive walls by two fictitious neutral walls at distance 
L'x^^ . We estimate L%^^ and p'^^^ in the following 
manner: for each configuration picked up for averages, 
we search the smallest and highest value of x and put 
Ll^^ = Xraax " Xmm + CT, and then p^ff = N/{L%f^Ly). 
A fixed value of L'!j^ can be obtained with different 
combinations of the parameters a and Ru,p in Uwp{x). 
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FIG. 2: Stripes morphology: (•) two stripes, (A) stripes perpen- 
dicular to the walls, (*) three stripes, as a function of the wall- 
particle potential range Rwp and of the effective wall separation 
L%^^. All simulations are performed with the nominal density 
p = 0.3 and T = 0.5, using N = 400 and = 32. The effec- 
tive densities, p*^^ ^ , range from 0.37 to 0.41. 



If Rwp{x) < Lx/4:, L%ff can be estimated as the length 
at which the foUowing relation holds: 



L^J^ /2) + [/* = T, 



(5) 



U* being the bulk energy per particle of p'^^f. The dis- 
crepancy between the values of L'^^^ obtained through 
eq. [5]and via simulations is less than 5 — 6% for the cases 
shown in Fig. [21 

Now, let us discuss in parallel the patterns obtained 
with true neutral walls as a function of (Fig. [1]), 
and those corresponding to fictitious neutral walls with 
L%ff = (Fig.EJ). First of aU, with = 26 the poten- 
tial energy profile has a minimum corresponding to the 
formation of three stable parallel stripes; with L'^Jf = 26, 
instead, we can stabilize parallel or perpendicular stripes 
depending on Ri^p (Fig. [2]). Interestingly, such a change 
occurs for R^p ~ i?r, i-e. the interparticlc repulsion 
range. As decreases towards 22, the system with 
true neutral walls moves towards a frustrated configura- 
tion having both parallel and perpendicular stripes. As 
L%^^ decreases towards 24, the stripes switch from paral- 
lel to perpendicular for smaller values of Rwp] finally at 
Ij<^ff 23 the parallel orientation is once again preferred, 
but it shows significant changes in the stripes width and 
period, with respect to the bulk and even to the case with 
neutral walls. In particular, in the region C of Fig. [2] A 
is not affected by the repulsion while P can vary up to 
20% with respect to the bulk; in the region B, P is not 
modified while minor changes of A up to 8% are present; 
in the region C both P and A are modified up to 25-30%. 
We conclude that strong confinement not only affects the 
orientation of the stripes but also the basic parameters, 
period P and A, of the microphasc. In this respect the 
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FIG. 3: (color online) Specific heat at constant volume per particle 
vs temperature for bulk and for different wall separations Lx- Cv 
is obtained as energy derivative. Confined systems are subject to 
neutral hard walls. Snapshots of the fluid at different temperatures 
are plotted in the insets. Periodic boundary conditions are applied 
along the y direction (dashed lines in the insets). 



behaviour of our system is rather different from the one 
reported in [2^ . in which asymmetric block copolymers 
are confined in a cylindrical pore of diameter D. When 
cylinders form within the neutral pore, those are par- 
allel or perpendicular to the wall depending whether D 
is a multiple of the bulk period or not. When the con- 
fining pore wall becomes attractive or repulsive a wide 
number of different morphologies is obtained in 27[ , the 
period of which is the same as in bulk. In our case, 
instead, the morphology at intermediate T remains es- 
sentially striped, but the values of P and A show some 
variation depending on the particle-wall interaction. A 



possible origin of such different behaviour in the two sys- 
tems under confinement is that in case of ref. p7j the 
diameter of the cylinders is much smaller of the period of 
the array whereas in our system A and P are comparable 
We find that also the temperature has a noticeable ef- 
fect on the microphase morphology under confinement. 
In bulk the onset of microphascs is marked by a peak in 
the specific heat. In systems confined between neutral 
walls, such a peak becomes broader and reduced in am- 
plitude, untill for = 14 it is barely a plateau. At the 
same time a secondary peak develops at low temperatures 
(T 0.36), for very narrow systems {L^ = 14), with a 
droplet-stripe passage driven by temperature (snapshots 
at low T in Fig.[3|). The results for = 24 suggest that 
such secondary peak in the specific heat is still present 
but shifted at a temperature below the lowest T of our 
computations. Likely such a peak disappears when 
is sufficiently large, since the system resembles more and 
more the bulk one, where the stripes are the most stable 
phase, as discussed in [3l|. Increasing the wall separa- 
tion, the droplets persist at low T mainly in the central 
region, while two stripes form close to the walls. The 
droplets are ordered as it is proved by the presence of 
Bragg peaks at short wave vectors in the static structure 



4 



factor. Circular clusters are entirely replaced by stripes 
for T > 0.4. 

In the bulk, the period slightly changes with the temper- 
ature [2^. Therefore one might suspect that the droplet 
formation is due to a change of P with T, so that a given 
Lx is commensurate to P at a certain T while it is not 
at another one. This is not the case as shown by low 
T simulations with ranging from 13 to 19: droplets 
always form at low T. We are therefore confident that 
such droplet formation at low T is not due to frustration 
like that governing the stripe orientation between neu- 
tral walls. In conclusion, the droplets can occur under 
two conditions: i) at fixed low T, droplets are favoured 
by strong confinement; ii) at fixed Lx, droplets appear 
at low enough temperatures. The role of energy and en- 
tropy in the stripe orientation as well as in the droplet 
formation will be discussed in a more extensive paper. 

In summary we have analysed a 2D microseparated 
fluid under lateral confinement, at a density such that 
stripes are formed in the bulk. We find that, with neu- 
tral walls, stripes perpendicular to the walls are favoured 
for Lx ^ P + A, P bulk period and A stripe width. 
Unfrustrated parallel stripes are obtained for wall sepa- 
ration Lx ~ nP + A {n integer). Otherwise frustrated 
configurations, made up of parallel stripes close to the 



walls and perpendicular stripes in the central region are 
stabilized, in contrast with the bulk case. The effect of 
a repulsive wall cannot be represented by a simple re- 
duction of the space L!^/-^ available to the particles. In 
fact for a given i^^/-^ , different patterns can be stabilized 
depending on the wall-particle potential range. This is 
not possible with neutral walls, since a unique configu- 
ration is selected for each separation Lx- Moreover, in 
presence of repulsive walls, a striped phase can be stabi- 
lized, characterized by rather different values of P and A 
with respect to the bulk ones. Finally, the role of tem- 
perature is important and droplets are favoured at low T 
for a system under confinement. 

Our results show that strong confinement is a way to con- 
trol the microphase morphology also in unexpected ways. 
The easiness in switching the stripe orientation at fixed 
T, by tuning the wall-particle interaction, might be an 
interesting property to investigate in materials which are 
also optically active. Finally it is of general interest to 
test the stability of the pattern morphology on the sam- 
ple temperature as well as on the degree of confinement, 
in order to develop new materials. 
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